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Mira variables in the Galactic bulge with OGLE-II data 
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ABSTRACT 

We have extracted a total of 1968 Mira variables from the Optical Gravitational Lensing Ex- 
periment II (OGLE-II) data base in the Galactic bulge region. Among them, 1960 are asso- 
ciated with 2 Micron All-sky Survey (2MASS) sources, and 1541 are further identified with 
Midcourse Space Exploration (MSX) point sources. Their photometric properties are com- 
pared with those of Mira variables in the Large and Small Magellanic Clouds. We have found 
that mass-losing stars with circumstellar matter are reddened such that the colour dependence 
of the absorption coefficient is similar to that of interstellar matter. We also discuss the struc- 
ture of the bulge. The surface number density of the bulge Mira variables is well correlated 
with the 2.2-/im surface brightness obtained by the Cosmic Background Explorer (COBE) 
satellite. Using this relation, the total number of Mira variables in the bulge is estimated to be 
about 6 x 10 5 . The log P-K relation of the Mira variables gives their space distribution which 
supports the well-known asymmetry of the bar-like bulge. 

Key words: stars: AGB and post-AGB - stars: mass-loss - stars: variables: other - Galaxy: 
bulge - Galaxy: structure - infrared: stars. 



1 INTRODUCTION 

Mira variables are believed to be in the last stage of Asymptotic Gi- 
ant Branch (hereafter AGB) evolution. They have larger amplitudes 
(e.g. > 2.5 mag in V-band) than other red variables, i.e. semi- 
regular variables or irregular variables. As a result of the short life- 
times of Mira variables, one needs a rather large population of stars 
to find a sufficient number of Mira variables to perform a statisti- 
cal study on them. This means that a large area should be observed 
repeatedly for a long enough duration to cover several pulsational 
periods, usually a hard task to do. However, the situation has been 
changed in the last decade by the gravitational microlensing search 
projects: Optical Gravitational Lensing Experiment (OGLE; Udal- 
ski, Kubiak & Szy manski, 1997), MA ssive Compact Halo Objects 
project (MACHO; Alcock et al. 2000), Microlensing Observations 
in Astrophysics (MOA; Bond et al. 2001), and so on. Researchers 
investigating variable stars are now enjoying the wealth of the data 
provided by these projects. 

The region studied most extensively is apparently the Large 
Magellanic Cloud: the di scovery of fi ve parallel sequences in 
the log P-K diagram by IWooJ fcOOOh was followed by many 
papers (Noda et al . 2002, 2004; Kiss & Be dding, 2003 2004 ; 
iGroenewesenl 120041 : Ita et al. 2004a, 2004 b; Frase r et alj|2005h . 
Following the c onclusion of iGlass & Lloyd Evans 1 198 lj) and 
iFeast et al] ll989t) that the scatter around period-magnitude rela- 
tion of Mira variables is small for bolometric or near-infrared mag- 
nitudes, papers focused on infrared properties of Mira variables and 
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those of other kinds of red variables. However, the Galactic bulge 
region has not re ceived full attention except for the Baade win- 
dows ISchultheis & Glas J200l|: Glass and Schultheis 2002, 2003; 
Schultheis, Glass & Cioni, 2004). Wray, Eyer & Paczynski (2004) 
analyzed OGLE variables over larger area of the bulge, but they 
only discussed small-amplitude variables and not Mira variables. 
Wozniak, McCowan & Vestrand (2004) also made a detailed anal- 
ysis of microvariability of Mira variables in the OGLE data, but no 
effort was dedicated to searching for infrared counterparts, which 
is indispensable to study the period-magnitude relation. 

In this paper, we collected Mira variables in the OGLE 
data and obtained their counterparts in the 2 Micron A ll- sky 
Survey (2MASS) point source catalogue jCurtietaljEobl) and 
the Midcourse Sp ace Experiment (MSX) point source catalogue 
(Egan et al. 2003). After we present data analyses to obtain a cata- 
logue in Section|2| we discuss their photometric properties (Section 
|3J and their space distribution in the bulge (Section|4j- 



2 DATA 

2.1 Extraction of Mira Variables from OGLE-II Data 

The Optical Gravitational Lensing Experiment-II (OGLE-II) re- 
peatedly obser ved 49 fields betwee n I = —11° and +11° from 
1997 to 1999 IWozniak etai]l2002t) . Galactic coordinates for the 
fields are given in Table □ Light curves (typically with 200-300 
data points) for 221 801 objects in these fields are available on the 
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Table 1. The 49 OGLE-II bulge fields. £ and b, galactic coordinate; Afjyiira. 
the number of OGLE-II Mira variables obtained in this paper; c(r a and o!dcc , 
offset between the coordinates in OGLE and 2MASS; Ax, the derived ex- 
tinction value. A^Mira with the superscript * indicates a couple of objects 
have overlaps with other fields (see Tablel3land text). 



Fid. I b A^Mira dRA <foec A K 

(deg) (deg) (arcsec) (arcsec) (mag) 



and a maximum-to-minimum amplitude (Al) for each variable. 
Small 6 indicates regular variation and a well-characterized period. 
We put the following criteria for Mira variables: 

P ^ 100 days, (1) 
< 0.6, (2) 
A/ ^ 1, (3) 

and obtained 1 1 207 candidates. All the light curves were further 
examined by eye. We found many fake detections that show very 
similar light curves but are fainter than the corresponding real de- 
tection at neighbouring positions, and rejected them. We also re- 
jected some other candidates that showed double periodicity which 
is often seen among variab les on the D sequence (Mira variables 
are on the C sequence; Wood 2000). Finally, we obtained 1968 
Mira variables. Table|2|shows the first 10 lines of the resultant ta- 
ble. An entry with Fld= / and ID= N corresponds to OGLE data 
named 'buLsc/ JV.dat4'. Their properties of variation (P, 8, A/, 
and mean magnitude < I >) are accompanied by identification 
with 2MASS and MSX point source, which will be discussed in the 
following sections. The full version of the table will be electroni- 
cally available in the online version of this journal. The number of 
Mira variables for each field, iViviira, is listed in TableQ 

2.2 Cross-identification with 2MASS 

Near-infrared counterparts of the OGLE-II Mira variable s were 
search ed for in the 2MASS point source catalogue ICurti et all 
2003). In order to adjust the systematic differences in positions be- 
tween the two data sets, we adopted the different offsets in right 
ascension and declination for different fields (c?r a and du cc , given 
in Table 0. With these offsets, we found 1960 matches within a 
limit radius of 1 arcsec. Fig.Qshows the differences in the right as- 
cension and declination. It is easily seen that most of the OGLE-II 
Mira variables have their near-infrared counterparts within a radius 
of 0.5 arcsec. 

Light variability causes main photometric uncertainties under 
discussion in this paper. Although the amplitude of light variation 
decreases with the wavelength for Mira variables, so me of them 
still h ave an amplitude of up to 1 mag in the /^-band l Gla ss et alJ 
1995; Whitelock, Marang & Feast, 2000). Therefore, mean mag- 
nitudes are usually used to investigate the infrared properties of 
Mira variables. As the 2MASS catalogue gives us single-epoch 
magnitudes, unfortunately, we should be aware of the possible er- 
ror caused by the difference between the mean magnitude and the 
single-epoch one. This error is large in comparis on with the differ- 
ence between photometric systems. According to lCarpented 1200 lh . 
for example, the difference is less than 0.1 mag between K s in the 
2MASS system and K in the Las Camp anas Observatory (LCO) 
system, in which svstem llta et alJ J2004aj) derived the log P-K re- 
lation (used in eq.[6). Therefore, we ignore the difference between 
K and K s magnitudes among various systems in the following dis- 
cussions. 
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Internet. 1 After we excluded false data with error flags, we appl ied 
the phase dispersion minimization method iStellingw erllll978|) to 
all the light curves. We used the same software as llta et alH 2004a) 
did. It gives us a period (P), a regularity indicator of variation (9), 

1 OGLE homepage i http://www.astrouw.edu.pl/ %>7Eogle/) 



2.3 Cross-identification with MSX 

Mid-infrared counter parts were sea rched for in the MSX point 
source catalogue jEgan et all 12003). leading to 1541 matches 
within a limit radius of 10 arcsec. The positional differences are 
plotted in Fig. [2] Among the identifications with 2MASS and/or 
MSX point sources, 30 Mira variables were found to be doubly 
detected because they lie in the overlapping regions between the 
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Table 2. The first 10 lines in the result table of Mira variables. This is a sample of the full version, which will be available in the online version 
of this journal. In case of null identification with 2MASS or MSX point source, 99.99 and 0.0 are put for magnitudes and colours, respectively. 
Flag= 1 indicates the double detection in the neighbouring fields (otherwise 0; see Table. |3*1. 



Fid. ID 


2MASS 


MSX6C 


P 


6 


AI 


< / > 


J-K B 


H-K s 




[8] 


Flag 








days 




mag 


mag 


mag 


mag 


mag 


mag 




1 98 


18022945-3024145 
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Table 3. Double detections in the overlapping regions between the neigh- 
bouring fields. 



5 ra (") 



5 Dec (") 



Figure 1. Histogram of the positional differences between OGLE and 
2MASS coordinates for the cross-identified Mira variables. 




Figure 2. Histogram of the positional differences between OGLE and MSX 
coordinates for the cross-identified Mira variables. 



neighbouring OGLE-II fields, as listed in Table[3] We put the flag 
1 in the last column of Table |2| if a variable has a counterpart in 
another field (otherwise 0). 

The y4-band (8.28 /J.m with the width of 3.36 /im) has the high- 
est sensitivity among the six photometric bands of the MSX survey. 
The majority of identified objects have a quality flag Q — 4 (Excel- 
lent) and all have flags better than Q — 2 (Fair/Poor) in the A-band, 
while more than two thirds of objects have a flag of Q = 1 (Limit) 
or Q = (Not Detected) in other bands. Therefore, we decided to 
use only the A-band flux in our discussions. The y4-band magni- 
tudes (deno ted as [81) were o btained using the zero magnitude flux 
of 58.49 Jy (Egan et al. 2003). Names and [8] magnitudes for iden- 
tified sources are listed in Table|2| The 427 objects without the MSX 
counterpart have relatively short periods (P < 400 days) and blue 
colours [ (H — K)o < 1 ], i.e. fainter K magnitudes according to 
the period- K relation and little excesses in mid infrared. They are 
naturally expected to be faint in the A-band. 
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Entry 2. 



2MASS 



MSX6C 



Fid. 


ID 


Fid. 


ID 


1 


3024 


38 


3446 


2 


3933 


33 


3377 


3 


878 


37 


857 


3 


2286 
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375 


3 


2695 


37 


2835 
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4124 
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2209 
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4235 
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2509 
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4522 


37 


4656 


3 


5929 


37 


5926 
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6181 


4 


4722 


3 


7804 


4 


6265 
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360 


39 


315 


5 


84 


44 


805 
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367 


44 


1066 


5 


3681 


44 


4608 


5 


3916 


44 


5045 


5 


4498 


44 


5702 


5 


4968 


44 


6453 


8 


297 


9 


210 


18 


1821 


19 


1540 


18 


5324 


19 


4649 


20 


1 


34 


1945 


21 


2270 


30 


1627 


21 


2895 


30 


2129 


21 


5977 


30 


5077 


30 


1821 


31 


159 


30 


4589 


31 


1759 


31 


793 


32 


877 


40 


1674 


41 


957 


40 


3769 


41 


3089 



18020050- 
18050058- 
17530423- 
17540566- 
17530457- 
17540599- 
17540524- 
17530343- 
17530367- 



2947556 
2839528 
3018592 
3009027 
3005208 
2956549 
2955533 
2952379 
2944057 



17540688- 
17550808- 
17495238- 
17495216- 
17495113- 
17495111- 
17495312- 
17495414- 
18233244- 
18073300- 
18073420- 
17584836- 
18005408- 
18005341- 
18005369- 
18015154- 
18015346- 
18025457- 
17513443- 
17513529- 



-2934463 
-3009531 
■3023328 
■3020497 
2954571 
2952354 
-2948492 
2944533 
2206346 
2724048 
-2649523 
2920304 
2902532 
2858004 
2836083 
2902159 
2840498 
■2854244 
■3320272 
■3251234 



G001.0639-03.5153 
G002.3762-03.5326 
G359.6485-02.0981 
G359.9034-02.2048 
G359.8450-01.9840 
G000.0783-02.1041 

G000.0255-01.8727 
G000.1497-01.8011 
G000.2718-01.9909 
G000.3989-01.9203 
G000.0047-02.4062 
G359.2301-01.5432 



G359.6381-01.2952 
G359.6717-01.2748 
G359.7288-01.2483 
G359.7876-01.2184 

G003.7556-03.4096 
G004.2579-03.1377 



G001.5985- 
G001.6689- 
G001.9870- 
G001.7114- 
G002.0268- 
G001.9396- 



02.9364 
02.8940 
02.7149 
03.1135 
02.9441 
03.2480 



G357.2979-03.1164 



2.4 Interstellar extinction 

The OGLE-II fields were chosen in the relatively low extinction 
regions in the Galactic bulge. Still they are located at low galac- 
tic latitudes where the interstellar extinction is generally strong. In 
order to estimate the extinction for each OGLE-II field, we mea- 
sured how much the red giant branch shifts along a reddening vec- 
tor in the colour-magnitude diagram. Comparing the position of the 
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< 0.4 




Figure 3. Comparison of the extinction values between Aj bv lSumiH2004t) 
and Ak by us. 



branch for each field with that for the Baade window NGC 6522 
field, the positional difference was interpreted as being the result 
of the extinction difference between the two fields. The extinction 
value for the NGC 6522 field itself was taken as Eb-v = 0.5 from 
ISchultheis et alJ 120041) . The procedures were adopted on both the 
J — K versus K diagram and the H — K versus K diagram. Slopes 
of the reddening vectors were taken from R ieke & Lebfskvljl985h : 



A K /E( J -K) = 0.659, 
A K /E{H-K) = 1.78. 



(4) 
(5) 



The extinction values obtained from both diagrams agree very well 
within ±0.02 mag for every field. The results are given in TableQ 
As shown in Fig. [3] our result s are c onsistent with the Ai values 
for OGLE-II fields obtained bv lSumil f2004). Although a metallic- 
ity gradient within the bulge can cause shifts on our estimates due 
to the dependence of the colour of giant bra nch, such a gradient was 
not detected and is considered to be small llbata & Gilmorel 1995b 
Frogel, Tiede, & Kuchinski, 1999: lRamirez et alJl200ol) . A differ- 
ence of 0.2 dex in [Fe/H] results in a shift of 0.03 mag in Ak- No 
reddening correction was applied to mid-infrared data since they 
are small and uncertain. 



3 DISCUSSION: PHOTOMETRIC PROPERTIES 
3.1 Completeness of the sample 

The saturation limit of OGLE-II observation is about / = 11. 
As a result, the extinction-corrected Iq magnitudes are limited to 
Io > 10 in Fig.|4] Some bright variables could be missing from our 
sample. However, fewer than 10 per cent are brighter than 1 1 mag 
in /-band magnitude distribution of late-type M giants toward the 
Baade window obtained by Blanco, McCarthy & Blanco (1984). 
Their distribution shows a similar decrease toward brighter region 
(/ < 12) as is shown by our histogram in Fig.|5] Furthermore, kine- 
matic studies showed these bright Mira variables to have a small ve- 



>ly foreg: 
kfcl987t 



Sharpies, Walker & 



locity dispersion, and they are probably foreground objects in the 
outer bulge or disc l Feast & Whi teloc 
Cropper, 1990). 

On the other extreme, no variables fainter than Io = 19 were 
detected. As shown in Fig. |4] the fainter boundary for the short- 
period Mira variables is clearly delineated, but the long-period Mira 
variables seem to extend to below the detection limit. While the 
reddest sources in our sample show (// — K)o ~ 2, some OH/IR 
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Figure 4. Iq magnitude plotted against log P (top) and (H~K)q (bottom). 
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Figure 5. Histogram of mean /-band magnitudes for the OGLE-II Mira 
variables. 



stars near the Galactic Centre have as large (H — K)o as 3 mag 
(Wood, Habing & McGregor, 1998). Such red objects have also 
been found among the IRAS sources in the Galactic bulge between 
7° < \b\ < 8° (Whitelock, Feast & Catchpole, 1991). These ex- 
tremely red objects are undergoing heavy mass loss which leads to 
large (H — K)$. In such cases, /-band flux should be reduced to 
below the detection limit. Fig. |6| compares the period distribution 
of our sample with that for the IRAS sample by Whitelo ck et alJ 
While their sample lacks short-period, and probably blue, 
Mira variables, it contains more significant component at P > 
700 days than our sample. Probably, some long-period red vari- 
ables are missing from our sample. In this regard, it is interesting 
that the OGLE-II objects with [8] > 1 Jy (filled portion in Fig.|6| 
have a more similar period distribution to the IRAS sources. 
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Figure 6. Histogram of p eriods for variables foun d in this work and in the 
7° < \b\ < 8° region bv lWhitelock et all <199ll) . In the bottom diagram, 
variables with MSX counterparts are hatched, and among them ones with 
[8] > 1 Jy are filled. 
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Figure 7. Colour-magnitude diagram (top) and colour-colour diagram (bot- 
tom) for the Mira variables. 



3.2 Mass-losing O-rich star or C-rich star ? 

Schul theis et ail <2004 compared the red variables in three regions, 
namely the NGC 6522 field, the Large Magellanic Cloud (LMC), 
and the Small Magellanic Cloud (SMC). In the NGC 6522 field, 
they noted the absence of variables with (H — K)o > 1, which 
are rather common in both Magellanic Clouds (MCs). They as- 
signed those red objects to carbon-rich (C-rich) stars which are 
absent in the Galactic bulge. Fig. shows the colour-magnitude 
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Figure 8. Histogram of the colours (H — K)q of Mira variables for the 
bulge (this work) and the both Magellanic Clouds llta et all2004bl) . 
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Figure 9. Distributions on the colour (H — K)o versus the amplitude AI 
diagrams for the three regions. 



and the colour-colour diagrams of the Mira variables in the OGLE- 
II fields. In contrast to the NGC 6522 field, our sample shows a 
red tail extending to (H — K)o = 2. However, Fig. [8] shows the 
fraction of such red variables is so small that the difference is not 
significant. Fig.|8|also presents the colour distribution of Mi ra vari- 
ables i n both MCs, the data for which were collected from Ita^yd 



(2004b) using the same criteria, i.e. As Schulthe is et al 
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(2004) mentioned, the distribution is clearly different among three 
galaxies. This is caused by the different combinations of oxygen- 
rich (O-rich) stars and C-rich stars, as is discussed in the following. 

Fig. |9| shows the relationship between the colour (H — K)o 
and the amplitude AJ. Ilta et alJl2004bT) argued that the dichotomy 
for the LMC is caused by the difference between O-rich stars 
(bluer) and C-rich ones (redder). A distinction between O- and C- 
rich stars is not apparent in Fig.|S| but is clearly seen in Fig.|5|as two 
separate branches. While the distribution of O-rich Mira variables 
forms a rather sharp ridge at 0.3 < (H - K) < 0.4 for the LMC, 
the ridge is wider and stretches to a redder colour for the bulge. 
This can be explained by the fact that the giant branch star with 
the higher metal content has the redder colour as shown for Galac- 
tic glo bular clusters (Frogel, Cohen & Persson, 1983: lFerraro et all 
2000). On the other hand, the peak disappears for the SMC. This is 
probably because the lower metal content for the SMC reduces the 
fraction of O-rich stars. 

The dichotomy in Fig. [5] for the LMC disappears for the 
bulge. The O-rich branch dominates, and the redder part looks ex- 
tended from the branch. This is in accordance with the fact that 
so far no bright carbon sta r at the tip o f the AGB has been found 
within the Galactic bulge (Blanco & Terndrup 1989). The top of 
Fie. 1101 shows that Mira variables in the redder part have longer 
periods. The middle panel and the bottom one plot the colour 
Kq — [8] and the mid-infrared magnitude [8], respectively, against 
log P. These panels show Mira variables with periods longer than 
log P — 2.55 have [8]-band excesses indicating thick dust shells 
produced by mass-loss phenomena. It is interesting that the log P- 
[8] relati on clearly has a turn-off in contrast to the logP-[12] re- 
lation by Whitelock et al. ( 1991) with only long period variables 
(logP > 2.4). The slope turns up to a larger one (—13) from a 
value similar to that of the log P-K relation (—3.59). 
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Figure 10. Relationships between period and colour (top, (H — K)q; mid- 
dle, Ko — [8]) and period- [8] magnitude relation. Slopes of a solid line and 
a dashed line are —13 and —3.59 respectively. 



3.3 The log P-K relation 

Fie. llll shows the \ogP-Ko diagram for the OGLE-II Mira vari- 
ables. The solid li ne is a log P-K relation of the LMC Mira 
variables <Itaetalj l2004a) after the correction of difference in 
the distance modulus. Crosses representing Mira variables with 
(H — K)o ^ 1 occupy a fainter region because of the circumstellar 
extinction in the if -band. Fig. ll2l shows the difference between Kq 
and the log P-absolute K relation: 

M K = -3.59 log P + 1.30, (6) 

for each variable, against (H — K)o. For the both MCs, it is clearly 
shown that Mira variables with (H — K)o *5 1 deviate from the 
constant values. The circumstellar extinction is responsible for this 
A'-band darkening because the reddened variables have longer pe- 
riods as those in the bulge do (see Fig. llOI . 

The distribution for the bulge seems different from those for 
the MCs and the deviations seem to occur even for bluer Mira vari- 
ables. The distribution for the bulge has a large spread compared 
with the MCs. This is caused by the depth of the bulge along the line 
of sig ht. As stars in the bulge have a rather wide range of metallic- 
ities IZoccali et ai1l2003l and references therein), there could also 
be an effect from some kinds of metallicity dependence. Although 
llta et alj 1200 4a) claimed they found the metallicity dependence of 
the zero point be tween both MCs, it is still a subject of controversy 
(see Feast 2004, for a review). 

Solid lines in Fig. 1121 are least-square fitted lines for Mira 
variables with (H — K)o ^ 1 after one 2-a clipping. It is in- 
teresting that their slopes are similar to the reddening coefficient 
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Figure 11. The period-magnitude relation for the bulge Mira variables. 
Crosses for Mira variables with (H — K)o 1, and filled circles for the 
others. 



A K /E(H - K) (1.67 ± 0.20 for the bulge; 1.67 ± 0.07 for the 
LMC; 1.73 ± 0.17 for the SMC). This indicates circumstellar mat- 
ters around Mira variables have optical characteristics similar to 
those of interstellar matters in the near infrared range (JHK). 
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Figure 12. The deviation from log P-Mk relation plotted against the 
colour (if — K)o for the three regions. Solid lines are fitted for the red 
part with (if - K) > 1. 



4 DISCUSSION: THE STRUCTURE OF THE BULGE 
4.1 Surface number density 

In Fig. 1131 the number of Mira variables for each OGLE-II field 
with —4° < b < —3° is plotted against the galactic longitude 
(cross). The amount of Poisson noise along the y-axis is illustrated 
on the right-hand side of the figure. The distribution is rather flat 
between — 3° < I < 3° and drops to the half maximum at around 
\l\ = 5°. We also confirmed this feature in the distribution of red 
clump stars (filled circles in Fig. 1 13i using the result obtained by 
Sumi (2004). This is because the surface brightness of the bulge 
is a so-called 'boxy' -type and the brightness does no t vary along 
b = — 3° between the region (see, for example, fig. 1 bv lDwek et alj 
119951) . 

We compared the surface number density with the surface 
brightness obtained by the Diffuse Infrared Background Experi- 
ment (DIRBE) onboard the Cosmic Background Explorer satellite 
( COBE). We used the Zodi-subtracted mission average map which 
is available on the Web. 2 In order to extract the surface brightness 
of the bulge (excl uding the disc), we basically followed the method 
of lWeiland etai]<1994l) . First, we need to know the interstellar ex- 
tinction. Between the Ak obtained in Section |2~4| and the COBE 
colours at the OGLE-II fields, we found the relation 
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Figure 13. The galactic longitude versus the number of Mira variables 
(cross) and that of red clump stars (filled circle) for each OGLE-II field 
with -4° < b < -3° 



as shown in Fig.^| From this equation, we obtained (J — K)o = 
0.8 as the intrinsic colour of the bulge population. It is a reason- 
able colour for old and metal-rich stellar populati on. For example, a 
recent population synthesis model by Mou hcine & Lancor] 1200 3l) 
predicts a stellar population with Z = 0.008 (or [M/H]= -0.4) 
would have such a colour over a wide range of ages between 16 Gyr 
and 2 Gyr. We regard the above equation as being applicable to ar- 
eas other than those observed by the OGLE-II survey, and used this 
relation to estimate the interstellar extinction. 

After the correction of the interstellar extinction, we sub- 
tracted the disc contribution. The disc is described in an exponential 
form: 



J(J,&)dtoc = J(0,6)di» 



(8) 



A K = 0.73 x [-2.5 log(/i. 2S //2.2) + 0.14], 



http://lambda.gsfc.nasa.gov/product/cobe/ 



(7) 



where lo(b) is the scale-hight for each galactic latitude, fitting to 
the observed surface brightness between 10° < \l\ < 45°. With the 
obtained disc components subtracted, the surface brightness of the 
bulge component for each field is obtained by averaging the values 
within 0.5° from the central coordinate of the field. Fig. ll5l shows 
the relationship between the thus-obtained surface brightness and 
the number of Mira variables for each field. 

The specific frequency of Mira variables per unit apparent 
brightness is 0.059 Jy from the slope of the figure. Assuming 
a distance of 8 kpc for the bulge, the frequency per unit A'-band 
absolute brightness is 0.094 MIy -1 . For metal-rich globular clus- 
ters, Froael & White loc3il99l determined the specific frequency 
of long-period variables to be 0. 1 1 MIy _ 1 (see their eq. 2). Consid- 
ering that they didn't include variables with P < 190 d, whereas 
they account for about 10 per cent of our sample, the frequency for 
the bulge is 25 per cent smaller than that for the globular clusters. 
This can be because the high metal content in the bulge enhances 
the mass loss and abbrevia tes the lifetimes of Mira variables as 
iFrogel & Whitelockl fl998) suggested. The total number of Mira 
variables were estimated to be 6 x 10 5 by integrating the surface 
brightness distribution over the \l\ < 10° and \b\ < 10° region. 



4.2 Radial distribution 

The period-magnitude relation of Mira variables has been widely 
applied as a distance indicator. A recent example is Reikuba 
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Figure 14. The relationship between COBETJIRBE colours and the extinc- 
tion values obtained from 2MASS data (see Section l2"4l . 
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Figure 15. The relationship between the surface brightness of the bulge 
obtained from COfi£/DIRBE and the number of Mira variables for each 
OGLE-II field. 



(2004) where the relation is used to obtain the distance to the 
galaxy NGC 5128 (Cen A). Furthermore, Lah, Kiss & Bedding 

(2005) explored the three-dimensional structures of the MCs. 
We expect that the structure of the bulge can also be traced 
by Mira variables. It is considered to be shaped like a bar in- 
clined from the Sun-Galactic Centre line of sight by 15-40° af- 
ter many works with various meth ods (distribution of luminous 
late-type stars. [Nakada et alJll99ll gas dynanucsjBjnney^yil 



1991 



1991 



Nakai 1992; surface brightness distribution, Blitz & Soersel 



Dwek et al 



1995; distribution and kinematics of SiO masers, 



Deguchi et al. 2002; peak of luminosity function due to red clump 
stars. IStan ek et alJll994l; iNishi vama et alJl2005h . Using 104 Mira 
variableZ 1wtoelock"^Catchpoler^^9^) also showed the asym- 
metry of the bulge. 

In Fig. 1161 the distribution of OGLE-II Mira variables are pro- 
jected on the galactic plane. Distances to individual variables are 
obtained by adopting the relation {(>}. As the very red objects tend 
to deviate from this relation, we used only stars with (H — K)o ^ 1 
to avoid the circumstellar extinction. Unfortunately, we had to use 
single-epoch magnitudes from the 2MASS catalogue, instead of av- 
erage magnitudes which are usually used for the distance estima- 
tion. The uncertainty in our estimate from single-epoch magnitudes 
causes the distance uncertainty to rise by up to ±2 kpc in the ex- 
treme case. Still we notice that there are far less stars at nearer 




Figure 16. Distribution of Mira variables projected on to the galactic plane. 
The Sun locates at (0,0) and the Galactic Centre at around (X = 8, Y = 0). 
Solid lines correspond to I = ±12°. 



region to the Sun on the negative longitude side than on the posi- 
tive side. This asymmetry should reflect the elongated shape of the 
bulge. It is expected that more dense surveys, i.e. toward a num- 
ber of galactic longitudes, will reveal the detailed structure of the 
bulge. It is also interesting to investigate the kinematics of the Mira 
variables and to combine them with the space distribution. 



5 SUMMARY 

We collected 1968 Mira variables in the bulge from OGLE-II data 
base. 1960 objects were cross-identified in 2MASS point source 
catalogue, and 1541 objects in MSX point source catalogue. In- 
terstellar extinctions for the near-infrared magnitudes of 2MASS 
were corrected by our estim ates of t he exti nction values Ak , which 
are consistent with Aj by ISumil 120041) . They are also consis- 
tent with those obtained with the surface brightness colours from 
COBEfDIKBE data. Based on the catalogue of Mira variables, 
some photometric properties were discussed. Our results include: 
(1) the separation between O-rich Mira variables and C-rich ones in 
the colour-amplitude diagram, (2) the preference for periods longer 
than 350 d of mass-losing stars demonstrated in the period-colour 
diagram and the period- [8] diagram, and (3) that those mass-losing 
stars are affected by their own circumstellar matter that have similar 
reddening vectors to that of interstellar extinction. 

We also discussed the structure of the bulge. The number of 
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Mira variables is proportional to the C0B£7DIRBE surface bright- 
ness, for which interstellar extinction and disc contribution were 
corrected. Using the coefficient of proportionality, we estimated 
the total number of Mira variables in the bulge is about 6 x 10 5 . 
Our analysis provided a panoramic view of the bulge although it 
is a preliminary result due to the lack of repeated observations in 
the infrared. However, we clearly detected the asymmetry between 
the positive galactic-longitude side and the negative side, which 
is commonly considered as evidence of the bar-like shape of the 
bulge. 
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NOTE ADDED IN PROOF 

Croenewegen & Blommaert (2005) have also made an analysis of 
Mira variables in the OGLE-II data, and independently obtained 
some similar results to ours. In particular, they present estimations 
of the structure using simulations. 
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SUPPLEMENTARY MATERIAL 

The following supplementary material is available online. 

Table 2. The result table of Mira variables. In case of null 
identification with 2MASS or the MSX point source, 99.99 and 9.9 
are put for magnitudes and colours, respectively. Flag= 1 indicates 
double detection in the neighbouring fields (otherwise 0; see Table 
3). 



